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We study the internal structure of a static and spherically symmetric neutron star in the framework 
of an in-medium modified chiral soliton model. The Equations of State describing an infinite and 
asymmetric nuclear matter are obtained introducing the density dependent functions into the low 
energy free space Lagrangian of the model starting from the phenomenology of pionic atoms. The 
parametrizations of density dependent functions are related to the properties of isospin asymmetric 
nuclear systems at saturation density of symmetric nuclear matter po — 0.16 fm“®. Our results, 
corresponding to the compressibility of symmetric nuclear matter in the range 250 MeV < Kq < 
270 MeV and the slop parameter value of symmetry energy in the range 30 MeV < Ls < 50 MeV, 
are consistent with the results from other approaches and with the experimental indications. Using 
the modified Equations of State, near the saturation density of symmetric nuclear matter po, the 
extrapolations to the high density and highly isospin asymmetric regions have been performed. 
The calculations showed that the properties of ~ IAMq and ~ 2Mq neutron stars can be well 
reproduced in the framework of present approach. 
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I. INTRODUCTION 

Analysis of the neutron star structure is an interest¬ 
ing topic of the modern astrophysics. In particular, in 
the nuclear astrophysics the structure studies of neutron 
stars are related to the Equations of State (EOS) of nu¬ 
clear matter which describes the pressure density ver¬ 
sus energy density relation for a broad range of density 
values (for example, see recent review [l[ and references 
therein). The peculiarities of EOS are well know at the 
densities below the saturation density of symmetric nu¬ 
clear matter po while at the high density regions they are 
still remaining not clear. The high density behavior of 
EOS are poorly understood because of the difficulty of 
direct experimental accessibility in laboratories and be¬ 
cause of the absence of ab initio theoretical calculations. 
Therefore, from the experimental point of view, the neu¬ 
tron star studies may serve as a laboratory for under¬ 
standing the behavior of EOS at high densities. From 
the theoretical point of view, instead of ab initio calcula¬ 
tions one can start from the phenomenological framework 
taking into account the well known properties of EOS at 
the low density region and extrapolate to the high den¬ 
sity regions trying to describe the properties of matter 
under the extreme (high density and high temperature) 
conditions. 

In this context and if one able to formulate fur¬ 
ther in-medium modifications, a chiral soliton model of 
Skyrme, describing the single nucleon properties in free 
space [i y , or its variations including the explicit vec¬ 
tor mesonic degrees of freedom may serve as a start¬ 
ing point for the theoretical framework. The in-medium 
modifications may be expressed allowing the density de¬ 
pendencies of the constants entering into the initial free 
space Lagrangian. It is necessary to note, that in prin¬ 
ciple one should be able also to reproduce the medium 
dependencies in the effective Lagrangian starting from 


the first principles but it is not known yet the form of 
general low energy Lagrangian and the peculiarities of its 
ingredients. For this reason, in the present work we use 
an in-medium modified chiral soliton model Hi which 
may be considered as a truncated version of the general 
low energy Lagrangian. 

In Refs. [^, [^ the medium modifications were achieved 
putting a bit more phenomenology into the initial model, 
i.e. putting the density dependent functions into the free 
space Skyrme Lagran^n according to the pionic-atoms 
data at low energies [7| and properties of asymmetric 
nuclear matter at saturation density po- Although the 
in-medium modified Skyrme Lagrangians is assumed to 
be very truncated version of the possible general La¬ 
grangian, it must be applicable to the studies of nu¬ 
clear many-body problems in the spirit of chiral effec¬ 
tive Lagrangians. The pay for the truncation may be the 
possible deviations from the experimental observables in 
the sense of quantitative description. Nevertheless, the 
model has obvious virtues: i) it has the simplest La¬ 
grangian among the same class Lagrangians, and ii) it 
seams have all necessary ingredients for the qualitative 
description of the strong interaction physics. These ideas 
were the basic ruliM philosophy behind of the approach 
developed in Refs. [alSl] and we continue our model stud¬ 
ies in the present work. 

The model is phenomenological one and must pass as 
much as possible tests on its applicability to strong inter¬ 
acting systems comparing with other approaches and the 
experimental indications. The previous nuclear matter 
studies showed that the in-medium modified Skyrme 
term is responsible for preventing the collapse of nuclear 
matter to the singularity at high densities in analogy to 
the free space case where the Skyrme term is responsibly 
for the stabilization of finite size solitons. The modifi¬ 
cations showed that at some values of model parameters 
the properties of infinite and isospin asymmetric nuclear 
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matter can be reproduced well near the saturation point 
of symmetric nuclear matter po- 

It will be interesting to test the applicability of model 
to the strong interacting systems under the extreme con¬ 
ditions extrapolating the modified Equations of State to 
the high density regions. Therefore, in the present work 
we make further check of the basic philosophy consider¬ 
ing an application of the model for the studies of neutron 
stars structure. 


II. FORMALISM 
A. The Lagrangian 


Our starting point is the in-medium modified Skyrme- 
model Lagrangian described in Ref. @ 

r = (1) 

1^*2 = ^ arTr {doUdoU^) - ^ a,Tr {d,Ud,U^) , 

amTv {2-U-U+) , 

£e = m^aeSabsTr (Tall) Tr (ndoU^) , 

where Einstein’s summation convention is always as¬ 
sumed (if not specified otherwise). The chiral SU{2) ma¬ 
trix U = exp(2irQ7ra/En-) is defined in terms of the Carte¬ 
sian isospin-components of the pion field tt^ (a = 1, 2, 3). 
The density functionals entering into the Lagrangian (a,-, 
0‘s, Ct, (s, a-m and ae) represent the influence of sur¬ 
rounding environment to the single soliton properties. 

In the single skyrmonic sector, the input parameters of 
the model have the following values: F,r = 108.783MeV, 
e = 4.854 and = 134.976 MeV. They correctly re¬ 
produce the experimental values of the nucleon mass 
Tow = 938 MeV and A-isobara mass uia = 1232 MeV 
in free space, i.e. if ar = otg = ctm = Cr = Cs = 1 and 
ae = 0 . 

This simple Lagrangian describes the properties of nu¬ 
cleons in free space, the properties of nucleons in nuclear 
matter as well as the properties of isospin asymmetric 
nuclear matter starting at the same footing. The formal¬ 
ism for the classical solitonic solutions, the quantization 
method and the explanations of applications for sym¬ 
metric and asymmetric matter properties can be found 
in Refs. [^, and references therein. But for seflconsis- 
tency, in the next subsection III Bl we briefly outline how 
an infinite size nuclear systems can be described in the 
framework of the present approach and represent some 
final formulas. 


B. Nuclear matter 


In the thermodynamic limit at zero temperature, for 
a system of an infinite number of baryons uniformly dis¬ 
tributed in inhnite volume but keeping the density per 
unit volume constant, the binding energy per nucleon 
can be represented as 

£(A, (5) = £y(A)-f ey4(A, (5), (2) 


where £y and ea are known as volume and asymme¬ 
try energies, respectively. For the convenience, here 
we introduced the isoscalar A = p/po and isovector 
6 = 6p/p density parameters in terms of the isoscalar 

P — Pneutron 4” Pproton and isOVeCtOr 6p = Pneutron Pproton 

nuclear densities, and the normal nuclear matter density 
Po- 

In the framework of present approach, using the spher¬ 
ically symmetric approximation for a single soliton prop¬ 
erties via the hedgehog ansatz U = exp{iririF(r)/r} and 
considering an isospin asymmetric and infinite size nu¬ 
clear environment, one can get the symmetric and asym¬ 
metric parts of the binding energy per nucleon Q 

£y(A) = /ito(/ 2) + - TOAf, (3) 

oa{\6) = ^ ( 4 ) 

Here the functionals m, A and A 2 are defined as 


to(/2) 


A 2 

A 4 



00 

f (r,2 I sin F 

3e3E, J\^ 


sin^ F dx . 


(7) 


In the functionals above, using an infinite nuclear matter 
approximation, the initial medium functionals in Eq. o 
were rearranged by defining the new medium functions 
/i, 2 , 3 , 4 - The rearrangements are made in the following 
way 


1 -I- Cip — fi 
1 -I- C 2 P = /2 
1 + Csp = fs 

C 4 X 


l + CsA 


S =U = 


{asHs ’ 

0^7" 

- “e X-1 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 
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and the basic principle behind the linear density depen¬ 
dent parametrizations of functions fi was the simplicity 
of form. 

In the last four equations Ci (i = 1,5) are the model 
parameters and, therefore, the present approach is a sim¬ 
ple 5-parametric solitonic model of nuclear matter.^ As 
soon as we define 5-parameters at some specific density 
(e.g. the normal nuclear matter density, po) the Equa¬ 
tions of State for the symmetric and asymmetric nuclear 
matter are defined for any values of nuclear matter den¬ 
sity if one assumes that the extrapolations by means of 
the linear on density functions fi are valid.^ In the next 
subsection III Cl we concentrate our attention on the prop¬ 
erties of baryonic systems at high densities discussing the 
neutron stars. 


C. Neutron star 


As the first approximation in describing the neutron 
stars, one can consider a spherically symmetric and static 
mass distribution. Then any part of neutron star mass 
M (r) inside a sphere with radius r is given by the integral 

Ai{r) = Att f drr^f(r). (12) 

Jo 

Here £{r) is the mass-energy density distribution of the 
neutron star in the radial direction. Consequently, the 
total gravitational mass M of the neutron star with ra¬ 
dius R is defined by the condition 

M = M{r>R). (13) 


Further, in the spherically symmetric approximation, the 
pressure density change in radial direction inside the neu¬ 
tron star is given W the Tolman-Oppenheimer-Volkoff 
(TOV) equation [1,[^ 


dP{r) 

dr 


G£{r)M{r) 


1-b 


P{r) 

£{r) 


1 - 


1 + 


2GM(r) 

r 

47rr^P(r 

M(r) 


■ (14) 


Here G = 6.707 x 10~^^hc ^ is gravitational con¬ 

stant and P{r) is pressure density in radial direction. 
The boundary conditions for the functions entering into 
the equation are 


7W(0) = 0 and £:(0) 

— ^cent ■ (15) 

After solving TOV equation one can find the radius R of 
a star with central energy density ^cent- It is defined by 


^ We would like to note, although at present work we discuss 
the properties of an asymmetric and infinite nuclear matter the 
model can be applied for the studies of finite nuclei properties 
after fitting the density parameters Ci. 

^ Other details can be found in Refs. Hi. 


the pressure zero condition at the surface of the neutron 
star, P{r — R) — 0. 

To obtain the numerical solution for the profile of a 
star, one solves Eqs. m and (imi using the Equation of 
State 


P = P{£). (16) 

To find P{£) dependence in present approach, we note 
that the pressure and energy dependencies on the density 
parameter A for the neutron matter (<5=1) are given by 
equations 

P(A) = poA^^^, (17) 

£{\) = [e(A, 1)-|-TO7v]Apo- (18) 

where £(A, 1) is binding energy per nucleon in neutron 
matter. The system of parametric equations, Eq. (ED 
and Eq. dUD, gives the desired relation P{£) between 
pressure and energy densities in a neutron star. 


III. RESULTS AND DISCUSSIONS 


In general, the thermodynamically limiting binding en¬ 
ergy per nucleon given in Eq. m must be valid at any 
densities and at any given isospin parameter S. For ex¬ 
ample, it is valid also at the extreme conditions which 
are formed in interior regions of neutron stars. Those ex¬ 
treme conditions described by very high density {p sev¬ 
eral times of po) and highly isospin asymmetric {6 cx 1) 
form of nuclear matter. 

From other side, it is not clear the direct relation of 
the liquid-drop formula of Bethe and Weizsacker [13, [HI 


B{Z,N) = avA-asA^/^-ac-pTR-o-A^ ^ 


Al/3 


A 


(19) 


describing the binding energy of nucleus to the bind¬ 
ing energy of neutron star. The reason is not only due 
to the relativistic metric factors coming from the Ein¬ 
stein’s equations in the calculations of binding energies 
of neutron stars. The reason here is rather obvious, the 
semiempirical liquid-drop formula describes the binding 
energy per nucleon near the normal nuclear matter den¬ 
sities (which correspond to the density profiles of the ex¬ 
isting heavy nuclei) and its validity at high densities is 
not clear. 

Nevertheless, around the normal nuclear matter den¬ 
sity and for the small values of asymmetry parameter d, 
the thermodynamically limiting binding energy per nu¬ 
cleon given in Eq. @ is well related to the Bethe and 
Weizsacker’s formula. Because, from one side, in the 
limit of small S the binding energy formula in Eq. @ 
can be approximated as 


e(A,d) = ey(A)+es(A)(52 + 0(^4)^ (20) 
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where es(A) is called the symmetry energy. From other 
side, if one ignores the Coulomb and surface effects, the 
binding energy per nucleon defined from the liquid-drop 
formula will take the form 

A-^B{Z,N) Ki av-aAS"^ + .... (21) 

From the comparisons of the last two equations, it is seen 
that the density dependencies of volume and symmetry 
energies can be established well at the densities around po 
from the stability conditions and the density variations 
of resonating heavy nuclei near the ground state. 

The symmetry energy eg describes the energy increase 
in the system if the number of protons and neutrons be¬ 
comes not equal relatively to the case when the neutron 
and proton numbers are same. Consequently, the symme¬ 
try energy is important factor in describing the properties 
of neutron-reach stable nuclei existing in nature as well 
as the properties of exotic nuclei formed under extreme 
conditions where the neutron-to-proton number N/Z is 
much smaller or much larger than one Although 

the definition of symmetry energy in Eq. dsni) is model 
independent its density dependence is clearly model de¬ 
pendent. Therefore, during the fitting to density region 
near the normal nuclear matter density po the different 
models are modified taking into account the properties 
of symmetry energy coming from the phenomenological 
observations. But the extrapolations of EOS to the high 
density (p >> po) and highly asymmetric {S ~ 1) re¬ 
gions remain not clear. In particular, this is due to the 
reason that in neutron stars the higher order terms in <5 
in Eq. (Uni) may be also important. Usually, it is assumed 
that 0(6'^) terms are small and the symmetry energy is 
mostly responsible in describing the properties of neutron 
stars when the neutron-to-proton number N/Z becomes 
infinite.^ In the present approach, the higher order terms 
0(6'^) come from the explicit isospin breaking symmetry 
in the mesonic sector and found to be negligible.^ Con¬ 
sequently, we ignore them in the present work. 

A. Nuclear matter and Symmetry energy 

Let us first discuss the properties of the symmetric nu¬ 
clear matter. The value of the first experimental param¬ 
eter in the liquid-drop formula Eq. (USD is well known, 
ay ^ 16MeV. It is also well established that the satu¬ 
ration density of symmetric nuclear matter, where the 


^ Although the neutron star is initially formed from the ordinary 
matter which has the finite neutron-to-proton number with well 
separated electrons, due to the gravitational collapse and due to 
further formed electron degenerate states, and due to the fol¬ 
lowing nucleon degenerate states at even higher densities the 
frequent collisions lead to the intensive nuclear reactions. As 
a result ”an effective” neutron-to-proton number N/Z becomes 
infinite. 

^ See the discussions in subsection V.C of Ref. [^. 
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FIG. 1: (Color online) The volume energy as a function of nor¬ 
malized density A = p/po- Solid, dashed and dotted curves 
correspond to the values of compressibility Kq — 230 MeV, 
Ko — 250 MeV and Ko = 270 MeV, respectively. Other pa¬ 
rameters are defined in Table H] 

pressure becomes zero Psym = 0 , is around 0.16fm“^. 
Further, the compressibility of symmetric nuclear mat¬ 
ter within the various approaches is found to be Kq ^ 
290±70MeV at the saturation density po [l3l - [l8l| . There¬ 
fore, these three quantities may be used to fix some of the 
model parameters (e.g. Ci, C 2 and C 3 ) if one expands 
the volume energy into the Taylor series 

ey(A) = ey(l)-I-^(A - I)^-I-... (22) 

around the saturation density, A = I. After fitting the 
values of parameters, the EOS of symmetric matter is 
defined for any values of density.^ The density depen¬ 
dencies of the binding energy per nucleon in symmetric 
matter are shown in Fig. [T] for the three possible values 
of the compressibility Kq. The corresponding parameters 
and the volume energy coefficients at saturation density 
Po are given in Table HI 

It is seen, that our results become consistent with the 
results from other approaches and the phenomenological 
indications at some values of density parameters Ci. For 
example, the density dependence of volume energy from 
the Set II is close to the APR (Akmal-Pandharipande- 
Ravenhall) predictions (see the model A18-|-(5u-|-UIX* in 
Ref. [ 1 ^) made on the basis of Argonne uis two-nucleon 
interactions . Similarly, the results from the Set III is 
very close to the results from the model A18-I-UIX M- 

The quantity in the last column, referred as the skew¬ 
ness parameter, is proportional to the third derivative of 


^ In this case, the minimization scheme is related to the single in¬ 
medium soliton, i.e. the soliton mass is minimized at the given 
values of density functions 
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TABLE I: Three sets of parameters which reproduce the sym¬ 
metric matter properties. The parameters C\, C 2 and Cs are 
chosen in such a way that at saturation point Psym(po) = 0 
the value of volume energy per nucleon equals to its experi¬ 
mental value, ev{po) = ~ —av, and the compressibility 

of nuclear matter Kq is reproduced in a given experimental 
range. 


Set 

Cl 

C 2 

C 3 

ev{po) 

[MeV] 

Ko 

[MeV] 

Q 

[MeV] 

I 

-0.285 

0.803 

1.753 

-16 

230 

-545 

II 

-0.273 

0.643 

1.858 

-16 

250 

-279 

III 

-0.333 

0.281 

3.090 

-16 

270 

-133 


the volume term at saturation density po and defined as 


Q = 27 


d^eviX) 

5A3 


A=1 


(23) 


Its values presented in Table U are outcome results from 
the present approach. There is a nonlinear correlation be¬ 
tween Q and compressibility Kq, Q increases if the value 
of Kq increases. Our predictions for Q is qualitatively 
similar to the results from the Hartree-Fock approach 
based on Skyrme interactions [2l|, and to the result 
from the MDI (isospin and momentum-dependent inter¬ 
action) model [22l| . Another example, the phenomenolog¬ 
ical momentum-independent model (MID) also predicts 
the similar results (23|. For comparison, in the present 
model one has QJKq ss —1.71 at Kq = 240 MeV while 
MID model gives the result Q/Kq Ki —1.6 at that value 
of Kq. 

Now let us discuss the properties of the asymmetric 
matter. As we said above, in the approximation that 
the higher order terms in 5 are small (see Eq. (l20l) l. the 
asymmetric matter properties are completely determined 
by the density dependence of the symmetry energy. The 
properties of the symmetry energy can be studied again 
by expansion into the Taylor series 


^5(A) = S5(l) + ^(A - 1) + §(A - 1)2 + ... (24) 


around the saturation density, A = 1. While the first co¬ 
efficient es(l) in Eq. (l24ll is known to be more or less well 
defined phenomenologically, es(l) ^ 29 to 34 MeV, the 
values of slop parameter Lg and compressibility of asym¬ 
metric matter Kg remain unclear. The reason is that 
EOS of asymmetric matter is highly sensitive to those 
parameters and the different models give the different 
predictions. For example, in relativistic mean field ap¬ 
proaches (see Ref. [ 2 ^, for the recent optimized versions) 
there are mainly two classes: i) small eg{pQ) ~ 30MeV 
and small Lg ~ 50 MeV (BSP, lUFSU*, lUFSU) and 
ii) large es(po) 37MeV and large Lg ^ llOMeV 
(Gl, G2, TMl*, NL3). The recent analyses of the data 
from heavy-ion collisions and neutron-skin experi¬ 
ments [ 2 ^ give the prediction Lg ~ 70 ± 20 MeV. 


During the modifications of EOS in the present work 
we choose the value of slop parameter Lg in the interval 
30—50 MeV which is in agreement with various empirical 
constraints [25l - [^ . 

The symmetry energies, calculated using two different 
sets of parameters, are shown in Fig. [2] and the parame- 
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FIG. 2: (Color online) The symmetry energy as a function of 
normalized density A = p/po- The solid and the dashed curves 
correspond to the sets ITa and Ill-f defined in Table UTl 

ters are determined in Table HIl Here we present only two 
representatives among the many sets producing the sym¬ 
metry energy parameters in commonly adopted range. 
Depending on the compressibility Kq value of the sym¬ 
metric nuclear matter we classify the results into two 
models. Model II with relatively smaller value of the 
compressibility Kq = 250 MeV (more soft EOS) and 
Model III with relatively bigger value of the compress¬ 
ibility Kq = 270 MeV (more stiff EOS). 

The results show that the symmetry energy is less sen¬ 
sitive to the different model parameters presented in Ta- 
bleini In the Eig.[2]we have shown the density dependen¬ 
cies of symmetry energy corresponding to two boundary 
regions referred as more soft (ITa) and more stiff (HTf) 
Equations of State. The other sets reproduce the density 
dependence of symmetry energy corresponding to some¬ 
where in between of these two boundary curves. The 
relatively more stiff EOS (HTf) in the present approach 
reproduce the density dependence of symmetry energy 
which is close to the APR predictions [1^, to the result 
from the MDI model with the parameter a; = 0 [1^ as 
well as to the results from MID model with the parameter 
2 /= -0.73 [ 13 . 

The quantities in the last five columns in Table HIl are 
outcome from our model calculations. In particular, the 
quantities Kj- and A"o ,2 are related to the compressibility 
of asymmetric matter and defined as 

Kr = Kg- 6Ls, 7fo.2 =Kr--^Lg. (25) 

Ao 
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TABLE II: The different sets of parameters which reproduce the asymmetric matter properties. The parameters Ci and Cs are 
chosen in such a way that, at saturation density pQ, the value of symmetry energy es and the coefficient of its first derivative 
Ls are reproduced in the commonly adopted range. Other parameters are defined in Table |T] (see the parameters values 
corresponding to the same Kq values given in this table). 


Set 

Ko 

[MeV] 

C 4 

Cs 

es(po) 

[MeV] 

Ls 

[MeV] 

Ks 

[MeV] 

Kr 

[MeV] 

Ko,2 

[MeV] 

£s(0.1fm ®) 
[MeV] 

£s(0.11fm ®) 
[MeV] 

Il-a 

250 

2.338 

0.878 

30 

30 

-209 

-299 

-265 

24.26 

25.54 

Il-b 

250 

1.984 

0.594 

30 

40 

-209 

-329 

-285 

23.11 

24.55 

II-c 

250 

1.723 

0.384 

30 

50 

-197 

-347 

-291 

22.06 

23.64 

Il-d 

250 

2.559 

0.946 

32 

30 

-222 

-312 

-278 

26.11 

27.44 

Il-e 

250 

2.183 

0.660 

32 

40 

-226 

-346 

-301 

24.93 

26.44 

Il-f 

250 

1.904 

0.448 

32 

50 

-217 

-367 

-311 

23.86 

25.50 

Ill-a 

270 

2.670 

1.498 

30 

30 

-169 

-259 

-245 

24.57 

25.76 

Ill-b 

270 

2.179 

1.024 

30 

40 

-177 

-297 

-278 

23.33 

24.71 

III-c 

270 

1.831 

0.701 

30 

50 

-172 

-322 

-298 

22.22 

23.75 

Ill-d 

270 

2.980 

1.622 

32 

30 

-178 

-268 

-254 

26.46 

27.69 

Ill-e 

270 

2.425 

1.133 

32 

40 

-189 

-309 

-290 

25.20 

26.62 

Ill-f 

270 

2.044 

0.798 

32 

50 

-188 

-338 

-313 

24.05 

25.64 


They describe the correlations between the symmetry en¬ 
ergy coefficients and important for estimating the shift in 
compressibility value in asymmetric matter 

if(po, 5)=K^ + Ko,2S^ + 0(6'^). (26) 

The condition for the lowering of saturation density value 
in asymmetric matter leads to the constraint < 0. 
This constraint is fulfilled in the present work. 

The calculated values of Kr and Kq 2 are consistent 
with the results from other approaches. For example, the 
phenomenological momentum-independent model pre¬ 
dicts the range for the values of A' 0 . 2 : —477 MeV < 
Ko ,2 < —241 MeV [ 2 ^. Our predictions, made using the 
different sets, are also belong to that range (see Table ITll). 

It is interesting also to compare the low density be¬ 
havior of the symmetry energy in the present model with 
other model predictions. For example, an analysis of 
the giant dipole resonance (GDR) of ^°®Pb using a se¬ 
ries of microscopic Hartree-Fock plus Random phase ap¬ 
proximation calculations predicts the following values of 
symmetry energy at subnormal nuclear matter density: 
20MeV < es{p = 0.1 fm’^) < 25.4MeV (H. The recent 
analysis of the properties of double magic nuclei [ 2 ^ puts 
more constraints, es{p = 0.1 fm“^) = 25.4±0.8MeV. An 
analysis of data on the neutron skin thickness of Sn iso¬ 
topes and the isotope binding energy difference for heavy 
nuclei at slightly higher density valu e g ave the result 
es{p = 0.11 fm~^) = 26.65 ±0.2MeV [S^. For compar¬ 
ison, our results are presented in the last two columns 
of Table EH One can see that they are mainly consistent 
with the results in Refs. [2§ - l^ . 

The discussions of density dependencies of the pressure 
in symmetric and asymmetric matter within the present 
approach can be found in Ref. [^. Here we only would 
like to note, that the modifications are also consistent 


with the results from other approaches. Moreover, the 
density dependencies of pressure in symmetric and asym¬ 
metric matter are weakly sensitive to the change of model 
parameters and the reproduced values correspond to the 
allowed region deduced from the experimental flow data 
and simulations studies by Danielewicz et al. [3l|. 

In summary, the properties of nuclear matter at satura¬ 
tion density po and the extrapolations to the lower than 
po regions are found to be satisfactory and consistent 
with the results from other approaches. Therefore, we 
now concentrate on the extrapolations to the high den¬ 
sity regions and discuss the properties of neutron stars. 


B. Properties of neutron stars 

Let us start from the mass-radius relations of neutron 
stars. There are dramatic differences in predictions from 
the different approaches as concerned the mass-radius re¬ 
lation of neutron stars. But most of the models with 
Equations of State of normal (nonstrange) nuclear mat¬ 
ter predict existence of the regions around the one solar 
mass where the radius of neutron star is independent 
from the mass. 

In the Fig. [3] we present our results corresponding to 
the two possible values of the compressibility of sym¬ 
metric matter Kq = 250 MeV in the left panel and 
Kq = 270 MeV in the right panel, respectively. One can 
see that the mass independent regions exist also in the 
framework of the present approach. Another the common 
for many models peculiarity, the existence of the maxi¬ 
mal mass of the neutron star, is also reproduced well in 
this model. As we mentioned above, the results are sen¬ 
sitive to the compressibility value, if the compressibility 
value decreases it leads to decreased values of the max- 
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FIG. 3: (Color online) The mass-radius relations of neutron stars at the values of compressibility Kq = 250 MeV (left panel, Set 
II) and Kq = 270 MeV (right panel, Set III). The value of symmetry energy at saturation density po is chosen as £s(l) = 32 MeV. 
The solid, dashed and dotted curves correspond to the symmetry energy slop parameter L values 30, 40 and 50 MeV. They 
represented by sets d,e and f in the Models II and III, respectively (see Table ED). 


imal mass and radius of neutron star leading in such a 
way to more compact neutron stars. From other side, 
the results are not sensitive to the change in the value of 
es(l) and, therefore, we presented only the results cor¬ 
responding to the value of £s(l) = 32MeV. In general, 
our results are in qualitative agreement with the results 
from other approaches. 

It is also interesting to compare our results in quan¬ 
titative level too. Measurements of the thermal spectra 
from the quiescent low-mass X-ray binaries inside globu¬ 
lar clusters gave the possibility to fit the data sets with a 
neutron star radius i?NS = 9-ill's ^1 90% confidence 
level [ 3 ^. Determinations of the mass-radius relation, 
based on recent observations of both transiently accret¬ 
ing and bursting sources, gave the radius range between 
10.4 and 12.9 km for 1.4 solar mass neutron stars [s^. 

In the present approach some class of parameters 
(Model III) defined in Table ED give very good agree¬ 
ment with the above mentioned estimations in quanti¬ 
tative level. The properties of neutron stars reproduced 
using some subclasses of Model III is presented in Ta¬ 
ble EID One can see that, our results corresponding to 
more than 2Mq as well as 1.4Mq neutron stars are 
very similar to the estimations from the Refs. [H, 

For comparison, in Table ITIIl we present two of the pos- 
sible neutron stars parametrizations from the Ref. [2l| as 
the representatives of the Skyrme effective forces used in 
the density functional approach. One can see that our re¬ 
sults in qualitative agreement with the results from the 
Ref. [llj. It is also necessary to note, that the compress¬ 
ibility value used in Ref. is smaller Kq = 230 MeV 
in comparison with the compressibility value in Model 
III Kq = 270 MeV. Therefore, we have more stiff EOS 
leading to slightly higher maximal mass and larger ra¬ 
dius neutron stars. But all of the models presented in 


Table EID give more or less similar results as concerned 
1.4Mq mass neutron stars. 

The central number densities of maximal mass neutron 
stars in the present approach are around 6.5po and the 
corresponding central mass-energy densities are around 
2.44 X lO^^gr/cm^ (« 1300MeVfm“^). Our results are 
close to the results from nonrelativistic potential model 
approaches discussed in Ref. [s^. 

It is interesting also to compare the total baryon num¬ 
ber of neutron star in the present approach with the re¬ 
sults from other approaches. General Relativistic for¬ 
mula for the total baryon number of the neutron star is 
given by the integral 


. . 2G7W(r) 

A = An J drr^p(r) - 


- 1/2 


(27) 


One can find the radial dependence in the number density 
p{r) from the relation P = P{p/po) after finding the 
radial dependence of the pressure P — P{r). Our results 
are in qualitative agreement with results from Ref. [2l| 
(see Table EID- 

After calculations of the total baryon number, one can 
also estimate the binding energy Ef, of the neutron star. 
Due to the decrease of the gravitational mass Ei, is de¬ 
fined by the formula 


Eh = Attin — M , (28) 

where we used the mass of nucleon ttin in free space.® 
One can see that the calculated binding energies corre¬ 
sponding to 1.4 solar mass neutron stars are consistent 


Note, that the authors of the Ref. [2l| used 1/56 part of the ®®Fe 
atom mass in calculating the binding energy of the neutron star. 
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TABLE III: Properties of the neutron stars from the different sets of parameters (see Tables HI and [III for the values of parameters): 
ric is central number density, pc is central energy-mass density, R is radius of the neutron star, Mmax is possible maximal mass, 
A is number of baryons in the star. Eh is binding energy of the star. In the left panel we represent the neutron star properties 
corresponding to the maximal mass Mmax and in right panel approximately 1.4 solar mass neutron star properties. The last 
two lines are results from the Ref. [^. 


Set 

ric 

pc 

[10^®gr/cm®] 

R 

[km] 

Mmax 

[Mq] 

A 

[10®^ 

Eh 

[10®® erg] 

Uc 

[fm-®] 

pc 

[10®®gr/cm®] 

R 

[km] 

M 

[Me] 

A 

[10®^ 

Eh 

[10®® erg] 

Ill-a 

1.046 

2.445 

10.498 

2.226 

3.227 

8.721 

0.479 

0.861 

11.587 

1.402 

1.898 

3.503 

Ill-b 

1.045 

2.444 

10.547 

2.223 

3.216 

8.557 

0.471 

0.861 

11.772 

1.402 

1.895 

3.453 

III-c 

1.037 

2.424 

10.616 

2.221 

3.200 

8.397 

0.460 

0.832 

11.953 

1.402 

1.887 

3.339 

Ill-d 

1.047 

2.452 

10.494 

2.221 

3.213 

8.598 

0.481 

0.867 

11.619 

1.402 

1.893 

3.422 

Ill-e 

1.044 

2.440 

10.554 

2.218 

3.203 

8.495 

0.473 

0.858 

11.809 

1.403 

1.890 

3.384 

Ill-f 

1.040 

2.433 

10.609 

2.216 

3.189 

8.311 

0.464 

0.842 

11.992 

1.403 

1.887 

3.334 

SLy230a 

1.15 

2.69 

10.25 

2.10 

2.99 

7.07 

0.508 

0.925 

11.8 

1.4 

1.85 

2.60 

SLy230b 

1.21 

2.85 

9.99 

2.05 

2.91 

6.79 

0.538 

0.985 

11.7 

1.4 

1.85 

2.61 


with the estimations made on the analysis of detected 
neutrinos from SN1987A: = 3.8 ± 1.2 x 10®^ erg [35j |. 

In addition we present the binding energy per unit 
mass versus GM/Rc^ relation in Fig. U] for some rep¬ 
resentatives of neutron star models given in Table IIIIl 
For comparison we represent also the outcome from an 



GM/m 

FIG. 4: (Color online) Binding energy per unit mass Et/M of 
the neutron star as a function of GM / Rc^. The solid, dashed 
and dotted curves correspond to the sets IITd, Ill-e and Ill-f 
defined in Table im Thin and black solid curve represents an 
approximate relation Eq. (PH]) presented in Ref. [34| . 

approximate formula 

Eh 0.6GM O.bGMY^ 

M ~ Rc^ J ^ ^ 

suggested by Lattimer and Prakash [s^ . It is seen that in 
the present approach the small values of slop parameter 
Ls ^ 30MeV of symmetry energy give the results close 
to the parametrization in Eq. (Ei). 


IV. SUMMARY 

In summary, we discussed the application of the in¬ 
medium modified chiral soliton model to the studies of 
asymmetric matter properties and neutron star struc¬ 
ture. The symmetric and asymmetric matter equations 
of state where reproduced by very simple 5-parametric 
density approach to the single nucleon properties in the 
nuclear environment introducing the isospin breaking ef¬ 
fects in the mesonic sector of the model. After reproduc¬ 
ing the asymmetric matter properties near the satura¬ 
tion density of symmetric matter pQ we extrapolated the 
Equations of State to the high density and highly isospin 
asymmetric regions. Our primary goal was a crude qual¬ 
itative analysis of neutron star properties. Nevertheless, 
at some given set of parameter values, our results are 
very close to the predictions from analysis of the data 
compiled during the observation of neutron stars [12,1121 ■ 
In particular, the calculations showed that the properties 
of ~ 1.4M0 and ^ 2Mq neutron stars can be well repro¬ 
duced within the present approach. 

As an outlook for further studies we note that the in¬ 
medium chiral soliton model presented here can easily be 
extended to the studies of finite nuclei properties. This 
task can be realized in the local density approach for the 
density of environment surrounding the soliton under the 
consideration. Those studies are under the way. 
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